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\^ ' Abstract. We present the results of spectrophotometry and V, R, I, Ha CCD photometry of the blue compact dwarf 

(BCD) galaxy SBS 0940+544. Broad-band images taken with the 2.1m KPNO* and 1.23m Calar Alto** telescopes 
reveal a compact high-surface-brightness H II region with ongoing star formation, located at the northwestern tip 
ON ' of the elongated low surface brightness (LSB) main body of the BCD. Very faint, patchy emission along the main 

body is seen in the Ha image. High signal-to-noise 4.5m Multiple Mirror Telescope (MMT)*** and 10m Keck 
II telescope^ long-slit spectroscopy of SBS 0940+544 is used to derive element abundances of the ionized gas in 
the brightest H II region and to study the stellar population in the host galaxy. The oxygen abundance in the 
brightest region with strong emission lines is 12 + log(0/H) = 7.46 - 7.50, or 1/29 - 1/26 solar, in agreement 
CIh' with earlier determinations and among the lowest for BCDs. H/3 and Ha emission lines and H5 and H7 absorption 

lines are detected in a large part of the main body. Three methods are used to put constraints on the age of 
the stellar population at different positions along the major axis. They are based on (a) the equivalent widths 
of the emission lines, (b) the equivalent widths of the absorption lines and (c) the spectral energy distributions 
(SED). Several scenarios of star formation have been considered. We find that models with single instantaneous 
bursts cannot reproduce the observed SEDs implying that star formation in the main body of SBS 0940+544 
was continuous. The observed properties in the main body can be reproduced by a continuous star formation 
process which started not earlier than 100 Myr ago, if a small extinction is assumed. However, the observations 
can be reproduced equally well by a stellar population forming continuously since 10 Gyr ago, if the star formation 
rate has increased during the last 100 Myr in the main body of SBS 0940+544 by at least a factor of five. We 
also constrain the age of the reddest southern region of the main body, where no absorption and emission lines 
are detected. On the assumption of zero extinction in this faint region, the observed spectrum can be fitted by 
a theoretical SED of a stellar population continuously formed with constant star formation rate between 100 
Myr and 10 Gyr ago. If, however, a small extinction of C(H/3) ~ 0.1 is present in this region then the observed 
spectrum can be fitted by a theoretical SED of a stellar population continuously formed with a constant star 
formation rate between 100 Myr and 1 Gyr ago. However, the poor signal-to-noise ratio of the spectrum and 
large photometric errors preclude reliable determination of the age of the southern region. In summary, we find 
no compelling evidence which favors either a young or an old age of SBS 0940+544. 
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1. Introduction 

Blue compact dwarf (BCD) galaxies are characterized by 
ongoing intense star formation (SF) activity as evidenced 
by strong nebular emission lines superposed on a blue con- 
tinuum (e.g. Sargent & Searle [1970 ; Lequeux et al. 1979: 
Kunth & Sargent [1983| Izotov, Thuan & Lipovetsky |l994fl . 
They are compact and metal-deficient objects with oxy- 
gen abundances ranging between 1/50 and 1/3 solar (e.g. 



et al. (1994), without an element abundance determina- 
tion. 

B and R surface photometry and colour profiles of SBS 
0940+544 have been derived by Doublier et al. (1997) (see 
also Doublier, Caulet & Comte ( 1999] )). 



Izotov & Thuan 1999) an d pos sess large am ounts of neu- 
tral g as (Thu an & Martin |l98l| ; van Zee et al. |l998| ; Thuan 
et al. |l999b|) . 

Most of the BCDs studied so far have revealed ex- 
tended low surface brightness red stellar sheets underly- 
ing the blue star-forming regions (Thuan 1983 ; Loose & 
Thuan 1986 ). Some of these stellar sheets have been re- 
solved with the Hubble Space Telescope (HST) into red 
giant stars (Schulte-Ladbeck, Crone & Hopp 1998; Crone 
et al. 2000[ ) . This indicates the presence of an evolved low- 
mass stellar populations in the underlying host galaxies. 

The tiny fraction of galaxies with extremely low oxy- 
gen abundance in the ionized gas (less than 1/20 (O/H)©) 
has been suggested by Izotov & Thuan ( |1999D to be 
young galaxies, based on chemical element abundance ar- 
guments. Nearly ten such galaxies with good abundance 
determinations are known to date. The outer parts of these 
galaxies are blue. The spectral energy distributions and 
colours of their outskirts are consistent with those of stel- 
lar populations with ages less than a few 1 00 Myr (see the 
examples of SBS 03 35-05 2 (Izotov et al. 
Izotov & Lipovetsky 1997 , Papade ros et a l. 
1415+437 (Thuan, Izotov & Foltz |l999a[ )) 



Correcting the 21 cm heliocentric velocity V = 1638 
km s -1 for Virgocentric infall with a velocity of 220 km 
s -1 and adopting a Hubble constant of 75 km s -1 Mpc -1 , 
Thuan et al. (1999b) derived a distance of D = 26.7 
Mpc for SBS 0940+544. At this distance, assumed here 
throughout, 1" corresponds to a projected linear scale 
of 129 pc. The total neutral hydrogen H I mass in SBS 
0940+544 is log(M(H i )/M ) = 8.71 (or M(H i) =5.1 x 
10 8 M ) (Thuan et al. |1999b[ ). 

In this paper we present new V, R, I and Ha photom- 
etry and very high signal-to-noise spectroscopic observa- 
tions of SBS 0940+544. We use those data to study the 
stellar populations in the BCD, put constraints on its age 
and derive elemental abundances in the ionized gas, in- 
cluding the helium abundance. In Sect. 2 we describe the 
observations and data reduction. The photometric prop- 
erties of SBS 0940+544 are discussed in Sect. 3. We derive 
elemental abundances in the brightest H n region in Sect. 
4. In Sect. 5 the stellar populations are discussed and con- 
straints on the age of SBS 0940+544 are considered. Sect. 
6 summarizes the results. 



2. Observations and data reduction 



1997H, Thuan, 2.1. Photometry 



1998) and SBS 



Detailed spectroscopic and photometric studies of ex- 
tremely metal-poor BCDs based on observations with 
large telescopes are useful to shed light on the origin of 
these galaxies, deduce their star formation history and 
constrain their ages. 

A good candidate for being a young galaxy is SBS 
0940+544. This galaxy was discovered in the course of the 
Second Byurakan Survey (SBS, Markarian, Lipovetsky & 
Stepanian 1983j ) and described as a compact extragalac- 
tic H II region without an underlying host galaxy, with 
strong emission lines, a weak continuum and an ultravio- 
let excess. 

A first abundance determination in the brightes t H H 
region of SBS 0940+544 was made by Izotov et al. ( |l99l|) 
who derived an oxygen abundance 12 + log(0/H) = 7.52 
± 0.12 and a helium mass fraction Y = 0.246 ± 0.030. 



Izotov, Thuan & Lipovetsky (|1997c|) , using higher S/N 
ratio observations, derived subsequently a lower value 12 
+ log(0/H) = 7.43 ± 0.01. Because of its very low metal 
content, SBS 0940+544 was used for the determination 
of the primordial helium abundan ce (Izo tov, Thuan & 
Lipovetsky |l994| ; Izotov & Thuan |l998b| ) and to study 
heavy element abundances (Thuan, Izotov & Lipovetsky 
1995; Izotov & Thuan 1999] ). The galaxy has been studied 



spectroscopically also by Augarde et al. (1994) and Comte 



Direct images of SBS 0940+544 were acquired in two ob- 
serving runs. Broad-band Johnson V and Cousins / im- 
ages and narrow-band images in the Ha line at A6606A 
with a full width at half maximum (FWHM) through 
a passband of 78A and in the adjacent continuum at 
A6520A through a passband with FWHM = 84A were ob- 
tained with the Kitt Peak 2.1m telescope on April 18 and 
April 22, 1999, respectively. The telescope was equipped 
with a lk Tektronix CCD detector operating at a gain of 
3 e - ADU -1 , giving an instrumental scale of 0'.'305 pixel -1 
and field of view of 5'. The total exposure of 40 and 60 min 
in V and /, respectively, was split up into four subexpo- 
sures, each being slightly offset from the others. Likewise, 
narrow-band exposures in the Ha line and the adjacent 
continuum bluewards of Ha were taken for 25 min each. 
The seeing during the broad- and narrow-band observa- 
tions was FWHM 1'.'2 and 2'.' 2, respectively. 

Another series of exposures with a total integration 
time of 2.5 hours in Cousins R was taken with the 1.23m 
telescope at Calar Alto in the period January 24 to 
February 22, 2000. The telescope was equipped with a 2k 
SITe detector at its Cassegrain focus giving a usable field 
of view of ~ 11' and an instrumental scale of 0V5 pixel -1 . 
The FWHM of point sources in the coadded R exposures 
is « T/5. 

Bias and flat-field frames were obtained during each 
night of both observing runs. Broad-band images were cal- 
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ibrated by observing different standard fields from Landolt 
(1992). Our calibration uncertainties are estimated to 



be well below 0.05 mag in all bands. Standard reduction 
steps, including bias subtraction, flat-field correction and 
absolute flux calibration were carried out using the IRAFQ 
and ESO MIDASpl. 



2.2. Spectroscopy 

High S/N long-slit spectroscopy of SBS 0940+544 was ob- 
tained with the 4.5m Multiple Mirror Telescope (MMT) 
and with the 10m Keck II telescope over a wavelength 
range 3600 - 7500A. 

The first set of observations was carried out with the 
MMT at two slit positions (see Fig. |l|b) . Observations were 
made with the blue channel of the MMT spectrograph op- 
erated with a 3072 x 1024 CCD detector. A l'/5 x 180" slit 
was used. The spatial scale along the slit was 0'.'3 pixel -1 , 
with a spectral resolution of - 7 A (FWHM). The final 
spatial scale is 0'.'6 pixel -1 after binning every two consec- 
utive rows into one. 

The first slit orientation with position angle P.A.=-41° 
(MMT #1, cf. Fig. lb), was centered on the brightest star- 
forming H II region (denoted as a in Fig. la) and oriented 
along the elongated main body of the galaxy close to the 
direction of the major axis. This 120-minute spectrum was 
obtained on December 13, 1996 at an airmass of 1.2. The 
total exposure time was broken into four 30 min exposures. 

The second slit orientation (MMT #2, cf. Fig. lb) is in 
the N-S direction with P. A. — 0°, centered on the second 
brightest region, denoted as c in Fig. [|a. These obser- 
vations were carried out on December 12, 1996 with an 
exposure time of 120 minutes at an airmass of 1.2. The 
total exposure time was also broken into four 30 min ex- 
posures. The seeing during both sets of MMT observations 
was 0'.'8. 

The Keck II observations were carried out on January 
9, 2000 with the low-resolution imaging spectrograph 
(LRIS) (Oke et al. [1995| ), using the 300 groove mm -1 
grating which provides a dispersion 2.52 A pixel -1 and 
a spectral resolution of about 8 A in first order. The slit 



was 1" x 180", centered on the brightest H n region and ori- 
ented with a position angle P. A. = -41°, the same as dur- 
ing the first set of MMT observations. No binning along 
the spatial axis has been done, yielding a spatial sampling 
of 0'.'2 pixel -1 . The total exposure time was 60 min, bro- 
ken into three 20 min exposures. All exposures were taken 
at an airmass of 1.2. The seeing was 0"9. 

Signal-to- noise ratios S/N ~ 50 and ~ 100 were 
reached in the continuum near H/3 for the brightest part of 



1 IRAF is the Image Reduction and Analysis Facility dis- 
tributed by the National Optical Astronomy Observatory, 
which is operated by the Association of Universities for 
Research in Astronomy (AURA) under cooperative agreement 
with the National Science Foundation (NSF). 

2 Munich Image Data Analysis System, provided by the 
European Southern Observatory (ESO). 



the galaxy during the MMT and Keck II observations. No 
correction for atmospheric dispersion was made because 
of the small airmass during all observations. This effect, 
although small, may introduce some uncertainties in the 
flux determination of the weaker lines. At an airmass of 
1.2 and position angle -41° the shift of the emitting region 
at the [O n] A3727 line relative to the same region at the 
[O in] A5007 line i n the direction perpendicular to the slit 
is 0"4 (Filippenko |l982| ). The shift at [O in] A4363 is much 
smaller, only O'.'l. In Sect. |we estimate the uncertainties 
introduced by this effect on the abundance determination. 
Spectra of a He-Ne-Ar comparison lamp were obtained 
after each exposure during the MMT observations, and 
spectra of a Hg-Ne-Ar comparison lamp were obtained af- 
ter each exposure during the Keck II observations. They 
were used for wavelength calibration. Two spectrophoto- 
metric standard stars, Feige 34 and HZ 44, were observed 
for flux calibration. The data reduction was made with 
the IRAF software package. The two-dimensional spec- 
tra were bias-subtracted and flat-field corrected. Cosmic- 
ray removal, wavelength calibration, night sky background 
subtraction, correction for atmospheric extinction and ab- 
solute flux calibration were then performed. 

One-dimensional spectra for the bright H II region a 
of SBS 0940+544 were extracted from the Keck II and 
MMT data within large apertures of 1" x 4" and 1'.'5 x 
3" respectively to minimize the effect of atmospheric dis- 
persion (Fig. ||). In addition we extracted spectra showing 
hydrogen Balmer absorption lines for four regions along 
the main body of the galaxy (i.e. along the slit oriented at 
PA. = -41°). The selected regions, denoted 1 to 4 (Keck 
II observations) and la to 4a (MMT) in Tables | and | 
were at adjacent positions along the major axis of the 
BCD, with the origin taken to be at the center of region 
a. Regions 1 and la coincide with region c in Fig. [lja. The 
spatial extents of the 4 regions along the slit were 3'.'0, 
2'.'4, 4"8 and 3'.'6 for the MMT observations and 3'.'0, 2"6, 
3"6 and 4"8 for the Keck II observations. Additionally, 
spectra with apertures 10"8 and 7'.'2 centered on region c 
were extracted from the two-dimensional MMT spectrum 
obtained at the position angle P.A. = 0° (regions 5 and 6 
in the Tables | and |) . 

3. Photometric analysis 

3.1. Morphology and colour distribution 

Star formation activity in SBS 0940+544 takes place pri- 
marily in the unresolved high-surface-brightness H II re- 
gion labeled a in Fig. [l]a. This region with an absolute V 
magnitude of -13.6 mag accounts for 80% of the total Ha 
emission and for nearly 1/4 of the total V light of the BCD 
within the 25 V mag arcsec -2 isophotc. A chain of faint 
point-like sources superposed on an asymmetric curved 
arm originating from region a is traceable to ~ 10" fur- 
ther to the northeast. This feature, labeled b in Fig. Qa, 
is not detected in Ha. Another feature is marginally visi- 
ble to the SW side of the elongated low-surface-brightness 
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Fig. 1. (a) V image of the cometary il BCD SBS 0940+544. The overlayed contours correspond to surface brightness 
levels between 23.5 mag arcsec -2 and 27 mag arcsec -2 in steps of 0.5 mag. The inset at the upper right shows a 
deconvolved 13" x 11" subframe at the northwestern tip of the galaxy. Adjacent to the brightest H n region a is 
a chain of other compact sources (labeled b) traceable to ~ 10" northeast of a. From the continuum-subtracted Ha 
map of SBS 0940+544 (right panel) and from spectroscopic study (Sect. 4), star-formation activity is found to be also 
present close to the optical maximum of the LSB body of the BCD (labeled c) . The average V — R and V — I colours 
of the main body increase from ~0.16 mag and ^0.43 mag in region c to ^0.34 mag and ^0.65 mag in region d. (b) 
Continuum-subtracted Ha map of SBS 0940+544 overlayed with V band contours between 21.5 mag arcsec -2 and 26 
mag arcsec -2 in steps of 0.5 mag. Roughly 80% of the Ha emission originates from the compact H II region a. The 
long-slit positions during the two observations with the Multiple Mirror Telescope (MMT #1 and MMT #2) and with 
the Keck II telescope (Keck #1) are indicated. 



(LSB) main body of the galaxy, at a surface brightness 
level <; 25.5 V mag arcsec -2 (region southwards of feature 
d in Fig. Ha). Judging from the continuum-subtracted Ha 
map (FigTfjp), star formation in SBS 0940+544 is not ex- 
clusively confined to region a but occurs also in the LSB 
host. A secondary Ha-emitting region is visible close to 
the geometrical center of the outer isophotes (region c in 
Fig. 0a). 

Evidence of recent star-formation activity along the 
major axis of the cometary LSB body of the BCD is also 
seen in the V — I and V — R colour maps (Fig. |J). The 
V — I and V — R colours in region c are respectively 
~0.43 mag and ~0.16 mag. The V — R colour index re- 
mains nearly constant at ~0.2 mag over a large part of the 
LSB body out to its southeastern extreme end (position 
+ 13",— 13" in Fig. ||b) and increases to ^0.34 mag at the 
faint (fiy ~ 24.0 — 24.5 mag arcsec -2 ) southern region d 
(position ~ +5", -14"). The V — I colour shows a simi- 
lar general, although less well defined trend, varying from 
~0.4-0.5 mag in the main body to ^0.65 mag in region d 
and getting bluer again (0.36 mag) close to the southeast- 



ern tip of the LSB host. The smooth colour change with 
decreasing intensity in the LSB body is accompanied by a 
clear isophotal twist between 23.5 and 27 V mag arcsec -2 , 
with the position angle P.A. changing from -40° to -10° 
and an average gradient of 8.5° mag -1 . Likewise, the el- 
lipticity 10 x (1 — ^) of the LSB host, a and b denoting 
respectively the semimajor and semiminor axis of the el- 
lipse fitted to a given isophote, decreases from 5.5 at 24 V 
mag arcsec -2 to 3.9 at 27 V mag arcsec -2 . 



Region a shows extraordinarily blue colours: V — I = 
-0.5 mag and V - R = -0.2 mag. Doublier et al. fll997|) 



have obtained B and R photometry of SBS 0940+544. 
They mistook region a for a Galactic foreground star. 
Spectra of the object show without doubt that it is an 
extremely compact H II region with very strong emission 
line s (Fig , g) at the redshift of SBS 0940+544. Izotov et 
al. (1991) first derived element abundances in this H n 
region, followed by Izotov et al. ( 1994| ) who named it SBS 
0940+544N. 
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Fig. 2. (a) V — I map of SBS 0940+544 displayed in the colour range -0.2 to +0.6 mag. The overlayed contours are 
from 21.0 to 25.5 V mag arcsec -2 in steps of 0.5 mag. (b) V — R map of the same region of the BCD displayed in the 
colour interval 0.1 to 0.4 mag. The overlayed contours are from 21.0 to 25.0 R mag arcsec -2 in steps of 0.5 mag. The 
bluest region in either map coincides with the northwestern H n region a, with an average V — I and V — R of -0.5 
mag and -0.2 mag, respectively. The angular resolution of each map is indicated by a circle, the diameter of which 
is equal to the FWHM of point sources in the corresponding smoothed input images used for its derivation (l'/8 and 
2'.'75 for the V — I and V — R map, respectively). 



3.2. Surface photometry and colours 



As remarked in Sect. 3.1 (see also Sect. [| and ||), signs of 
low-level ongoing or recent star formation are present all 
over the inner part of the main body of SBS 0940+544. In 
order to disentangle the stellar from the gaseous emission 
and better constrain the surface brightness and colour dis- 
tribution of its faint underlying LSB host, we performed 
surface photometry in all available broad-band filters in 
two steps. First we fitted the intensity distribution and 
subtracted both the dominant region a and all adjacent 
point sources (designated by b, cf. Fig. |l]a). We then de- 
rived surface brightness profiles (SBPs) for t he mai n body 
using methods described in Papaderos et al. ( |l996b[ ). SBPs 
computed in this way (Fig. |^a) trace the luminosity of 
the cometary LSB host of the BCD and the superim- 
posed star-forming zone extending from region c to the 
SE corner. The radial intensity distribution of the LSB 
host was approximated by adjusting different fitting laws 
to the outer part of the SBPs. We then computed SBPs for 
the whole BCD, and by subtracting from them the best 
model for the intensity distribution of the LSB host, we 
obtained the surface brightness distribution of the discrete 



and diffuse starburst sources located in the northwestern 
part and main body (cf. Fig. ||b). 



The SBPs of the main body (Fig. ||a) provide an up- 
per limit to the intensity distribution of the underlying 
unresolved older stellar population of the LSB host, since 
part of the light for R* 8" comes from region c and ad- 
jacent young sources. An exponential fitting law provides 
a reasonable approximation to the SBPs for /iy <; 24.5 
mag arcsec -2 , for radii < 10", however, it predicts a 
higher intensity than the one observed. This is also the 
case for the SBPs derived for the total emission of the 
BCD (Fig. |^b). This type of convex profile with an ex- 
ponential distribution in the outer parts, and leveling 
off in the inner part (within 1-3 scale lengths) is not 
rare among low-luminosity dwarf ellipticals (e.g. Binggcli 
dwarf irregulars (Patterson & Thuan 



fc Ca meron |1991|) 

1996| , Makarova |l999| , van Zee [2000| ) and blue low surface 



brightness galaxies (Ronnback & Bergvall 1994, Bergvall 
et al. 1999), and h as been reported i n a few BCDs (e.g. 
Vennik et al. 1996, Papaderos et al. 1996b, Telles et al. 



1997, Papaderos et al. 1999, Fricke et al. 2001). The ob 
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Fig. 3. (a) Surface brightness profiles (SBPs) of the main body of SBS 0940+544 after subtraction of the brightest H 
II region a and the adjacent point sources b (cf. Fig. |l|a). For the sake of clarity, the R and J SBPs have been shifted 
vertically by -1 and -2 mag respectively. The modeled surface brightness distributions of the underlying LSB host 
in the V band corresponding to a Sersic fitting law with an exponent r] — 2 (cf. Eq. and a modified exponential 
distribution with (b,q) = (3. 0,0. 9) (Eq. ||) are shown respectively by the thin/black and thick/grey curves, (b) SBPs of 
the total emission of SBS 0940+544. The meaning of the symbols is the same as in the left panel. Small open circles 
show the surface brightness distribution of the light in excess of the model for the LSB host (thick/grey curve) in the 
V band. The effective radius r e g of the V band SBP and the isophotal radii P25 and E25 of the starburst and LSB 
component at 25 V mag arcsec -2 are indicated. 



served profile shape suggests that a Sersic (|l990j) profile 
of the form 



I(R*) = J exp 



R* 



(1) 



with X] >1 provides a better approximation to the light of 
the LSB host than a pure exponential (77 = i). 

An alternative fitting formula which reproduces well 
an exponential distribution flattening inwards has been 
discussed by Papaderos et al. ( 1996b| ): 



R 



I(R*) = I exp ) [f - q exp(-P 3 (i?*))] , 

where Ps(R*) is 



bo.) \ a q 



(2) 



(3) 



Near the center, such an intensity distribution depends on 
the relative central intensity depression q = A// Jo, where 
Jo is the central intensity obtained by extrapolation of the 
outer exponential slope with the scale length a to R* = 0, 
and ba is the cutoff-radius where the central flattening 
occurs. 

Expressions (|]) and (||) were fitted to the SBPs of the 
main body for radii R* > 9", where the colours become 
nearly constant (Fig. |]). Fits with a Sersic profile give rj 
= 2.0 ± 0.1 and, as illustrated for the V profile in Fig. ||a, 



are nearly indistinguishable from those using Eq. g with 
b = 3.0 and q = 0.9. The V surface brightness distribution 
of the residual light in excess of the fit to the emission of 
the main body (Eq. ||) is shown by small circles in the 
right panel of Fig. ||. 

Table @ summarizes the derived photometric proper- 
ties of the LSB host and the star-forming component of 
SBS 0940+544. Cols. 2 & 3 give respectively the central 
surface brightness (ie,o an d exponential scale length a" 
of the LSB host as obtained from linear fits to the SBPs 
for R* > 9" and weighted by the photometric uncertainty 
of each point. The corresponding integrated magnitude 



2 lsb of a pure exponential distribution (ie,q—5 log(a")— 2 
is listed in col. 4. As discussed above, the surface bright- 
ness distribution of the LSB host becomes significantly 
shallower for R* < 9" than that predicted by inwards 
extrapolation of the exponential slope observed at the 
outskirts of SBS 0940+544. Therefore, for the case un- 
der study roLg B overestimates the apparent magnitude of 
the LSB host (w tolsb; col. 9). Cols. 5 through 9 list 
quantities obtained from profile decomposition, with the 
intensity distribution of the LSB host modeled by Eq. |[ 
Cols. 5 and 7 give respectively the radial extent P25 and 
E25 of the starburst and LSB components, determined at 
a surface brightness level of 25 mag arcsec -2 (cf. Fig. ||b). 
The apparent magnitudes of each component within P25 
and E25 are listed in cols. 6 and 8. Col. 9 gives the ap- 
parent magnitude of the LSB component in each band 
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Table 1. Structural properties of the starburst- and LSB component of SBS 0940+544 



Band 


fl>E,0 


a 


mtsB 


P25 


mp 25 


E25 


m E 2 5 


W-LSB 


msBP 


rratot 






mag arcsec -2 


arcsec 


mag 


kpc 


mag 


kpc 


mag 


mag 


mag 


mag 


kpc 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


V 


20.46+0.12 


2.51+0.05 


16.46 


0.82 


17.89 


1.36 


17.53 


17.34 


16.81+0.02 


16.83 


0.60,1.05 


R 


20.14+0.10 


2.52+0.05 


16.13 


0.76 


18.34 


1.46 


17.1 


16.97 


16.68+0.02 


16.67 


0.68,1.13 


I 


19.83+0.13 


2.50+0.06 


15.85 


0.83 


17.89 


1.54 


16.87 


16.76 


16.45+0.04 


16.48 


0.69,1.13 



within the photometric radius of 15'.'5, as derived from 
integration of the modeled distribution Eq. ^. The total 
magnitude of the BCD as inferred from SBP integration 
out to the same radius and by summing up the flux within 
a polygonal aperture is given in cols. 10 and 11, respec- 
tively. Col. 12 lists the effective radius r e g and the radius 
rgo which encircles 80% of the galaxy's total flux. 

The V — R and V — I colour profiles of the LSB host 
of SBS 0940+544 as derived by subtraction of the cor- 
responding SBPs (Fig. ||a) are shown in Fig. ||. Both 
colour profiles show very weak gradients for [iy <; 23 
mag arcsec -2 , and become nearly constant for R* > 9". 
Beyond that radius the colour profiles (Fig. 4) give aver- 
age colour indices V — R— 0.33 ± 0.04 mag and V — I 
= 0.58 ± 0.03 mag for the LSB component. These values 
are consistent with the integrated colours V — R = 0.37 
mag and V — I = 0.58 mag, determined from the appar- 
ent magnitudes of the modeled intensity distribution of the 
LSB host within R* = 15" (Table 0, col. 9). 

We have compared our results with those obtained for 
the B and R profiles by Doublier et al. ( [1997] , |1999[ ). The 
agreement is not satisfactory. The exponential scale length 
a of the LSB host derived here (Table [l]) is very similar 
in all bands (V R an d 7) being ~2'/5 (320 pc). Although 
Doublier et al. (|l997 ) also report an exponential intensity 
decrease in B and R for radii 3", they derive for the dis- 
tance assumed here a B and R exponential scale length 
of ~ 1 kpc, 3 times greater than our value. A striking dis- 
agreement is also evident from a comparison of the colour 
profile s derived here and those presented in Doublier et al . 
( |1997|) . The B-R colour profile of Doublier et al. ( |l997D 
shows a maximum of ~ 0.5 mag at R* ~ 0", then decreases 
monotonically with increasing radius to ^ 0.2 mag. Our 
colour profile shows the opposite trend: very blue colours 
at small radii and nearly constant, moderately red V — R 
and V - I colours in the outskirts of SBS 0940+544. 



4. Chemical abundances 

Earlier determinations of element abundances in SBS 0940 
+544 (Izotov et al 



1991| , |1994| ; Izotov & Thuan |l998b|) 



have revealed a metallicity ranging from Z@/31 to Zq/25, 
among the lowest for BCDs. Abundance determinations in 
this galaxy are therefore important for studying the origin 
of elements in low-metallicity environments and deriving 
the primordial helium abundance. We emphasize that re- 
peated observations of the same galaxy are very useful 
to estimate the robustness of the abundance determina- 
tions. This was one of the motivations for carrying out 



main body 




R* (arcsec) 

Fig. 4. V — R and V — I colour profiles of the main body 
of SBS 0940+544 derived by subtraction of the surface 
brightness profiles in Fig. |3|a. The radial range correspond- 
ing to the surface brightness interval between 22 and 23 
V mag arcsec -2 is indicated by vertical lines. The average 
colours of the bright H 11 region a (cf. Fig. 0a) ( V - R ~ 
-0.2 mag and V — I 55 -0.5 mag ) are indicated by hori- 
zontal lines. 



new MMT and Keck II observations of SBS 0940+544 
and analyzing the abundance pattern in its brightest H 
11 region. The spectra of region a used for the abundance 
determination are shown in Fig. ^. 

A two-zone photoionizcd H 11 region has been assumed 
for abundance determination; the electron temperature 
T e (0 in) for the high- ionization region has been derived 
from the observed flux ratio [O in] A 4 363/( A 4959+5007), 
using a five- level ato m mod el (Aller 1984 ) with atomic 
data from Mendoza ( 1983| ). The electron temperature 
T e (0 11) for the low-ionization region has been derived us- 
ing the empirical relation between T e (0 11) and T e (0 ill) 
from the H II region photoionization models by Stasinska 
( |1990| ). The [S 11] A6717/A6731 ratio was used to derive the 
electron number density A^ e (S 11). 

The spectra were corrected for interstellar extinction, 
where the extinction coefficient C(H/3) was derived from 
the hydrogen Balmer decrement using the equations given 
in Izotov et al. ( |1994|) and the theoretic al hy drogen emis- 
sion line flux ratios from Brocklehurst ( 1971 ). 

The emission line fluxes were measured using a 
Gaussian profile fitting. The errors of the line fluxes mea- 
surements include the errors in the fitting of profiles and 
those in the placement of the continuum. We also take into 
account the errors introduced by uncertainties in the spec- 
tral energy distributions of standard stars. Standard star 
flux deviations for both standard stars Feige 34 and HZ 44 
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Table 2. Observed (F(X)), extinction-corrected (/(A)) fluxes and equivalent widths (EW) of the emission lines in the 
spectra of the brightest H II region in SBS 0940+544 





Keck II 




MMT 




Ao(A) Ion 


F(X)/F(Rfi) /(A)//(H/3) 


EW (A) 


F(A)/F(H/3) I(A)/7(H/3) 


EW (A) 



3727 [O n] 

3750 H12 

3771 Hll 

3798 H10 

3820 [He i] 

3835 H9 

3868 [Ne III] 

3889 H8 + He i 

3968 [No in] + H7 

4026 He I 

4068 [S n] 

4101 HS 

4340 H7 

4363 [O in] 

4388 He 1 

4471 He 1 

4658 [Fe ill] 

4686 He 11 

4711 [Ar iv]+He 1 

4740 [Ar iv] 

4861 H/3 

4922 He 1 

4959 [O in] 

5007 [O in] 

5199 [N 1] 

5271 [Fe in] 

5876 He 1 

6300 [O 1] 

6312 [S in] 

6363 [O 1] 

6563 Ha 

6584 [N 11] 

6678 He 1 

6717 [S 11] 

6731 [S 11] 

7065 He 1 

7136 [Ar in] 

7320 [O 11] 

7330 [O 11] 

C (JS.fi) dex 
F(R/3) a 
EW(abs) A 



0.445 ±0.007 
0.013 ±0.001 
0.023 ±0.001 
0.042 ±0.001 
0.017 ±0.001 
0.058 ±0.001 
0.340 ±0.005 
0.184 ±0.003 
0.259 ±0.004 
0.014 ±0.001 
0.009 ±0.001 
0.244 ±0.004 
0.461 ±0.007 
0.125 ±0.002 
0.003 ±0.001 
0.036 ±0.001 
0.005 ±0.001 
0.006 ±0.001 
0.022 ±0.001 
0.013 ±0.001 
1.000 ±0.015 
0.010 ±0.001 
1.415 ±0.021 
4.251 ±0.061 
0.003 ±0.001 
0.004 ±0.001 
0.102 ±0.002 
0.011 ±0.001 
0.012 ±0.001 
0.004 ±0.001 
2.946 ±0.043 
0.013 ±0.001 
0.029 ±0.001 
0.038 ±0.001 
0.030 ±0.001 
0.026 ±0.001 
0.028 ±0.001 
0.007 ±0.001 
0.007 ±0.001 



0.471 ±0.008 
0.022 ±0.002 
0.032 ±0.002 
0.052 ±0.002 
0.018 ±0.001 
0.069 ±0.002 
0.356 ±0.006 
0.201 ±0.004 
0.278 ±0.005 
0.014 ±0.001 
0.009 ±0.001 
0.259 ±0.004 
0.477 ±0.007 
0.128 ±0.002 
0.003 ±0.001 
0.037 ±0.001 
0.005 ±0.001 
0.006 ±0.001 
0.022 ±0.001 
0.013 ±0.001 
1.000 ±0.015 
0.010 ±0.001 
1.402 ±0.021 
4.204 ±0.061 
0.003 ±0.001 
0.004 ±0.001 
0.097 ±0.002 
0.010 ±0.001 
0.011 ±0.001 
0.004 ±0.001 
2.750 ±0.043 
0.012 ±0.001 
0.027 ±0.001 
0.035 ±0.001 
0.028 ±0.001 
0.024 ±0.001 
0.026 ±0.001 
0.007 ±0.001 
0.006 ±0.001 

0.085 ±0.019 
1.21 ±0.01 
1.1 ±0.1 



62.0 ±0.3 
1.9 ±0.2 

3.3 ±0.2 
6.1 ±0.2 

2.5 ±0.3 

8.6 ±0.2 

49.6 ±0.3 
26.9 ±0.3 
38.3 ±0.3 

2.1 ±0.2 

1.4 ±0.2 
39.3 ±0.3 
84.9 ±0.4 
23.2 ±0.3 

0.6 ±0.2 
7.4 ±0.2 

1.2 ±0.3 

1.3 ±0.3 
5.2 ±0.3 

3.0 ±0.2 
241.1 ±0.8 

2.4 ±0.3 
348.8 ±1.0 

1052.0 ±1.7 
0.9 ±0.2 

1.2 ±0.3 
39.0 ±0.5 

4.9 ±0.4 

5.3 ±0.3 
1.9 ±0.3 

1409.0 ±2.5 
3.9 ±0.3 
14.9 ±0.5 

19.2 ±0.4 

15.3 ±0.4 

14.7 ±0.4 
16.5 ±0.5 

4.6 ±0.6 

4.1 ±0.4 



0.585 ±0.008 
0.022 ±0.003 
0.044 ±0.003 
0.059 ±0.004 



0.587 ±0.008 
0.023 ±0.004 
0.045 ±0.005 
0.060 ±0.005 



0.064 
0.373 
0.204 
0.288 
0.019 
0.007 
0.267 
0.487 
0.142 



±0.004 
±0.006 
±0.005 
±0.005 
±0.004 
±0.003 
±0.005 
±0.007 
±0.004 



0.065 
0.374 
0.205 
0.290 
0.019 
0.007 
0.269 
0.488 
0.142 



±0.005 
±0.006 
±0.006 
±0.006 
±0.004 
±0.003 
±0.006 
±0.007 
±0.004 



0.039 ±0.003 0.039 ±0.003 
0.003 ±0.002 0.003 ±0.002 



62.7 ±0.7 
2.5 ±0.3 
4.9 ±0.4 
6.7 ±0.5 

7.4 ±0.4 

44.0 ±0.7 

25.1 ±0.6 

39.8 ±0.7 
2.7 ±0.5 
1.0 ±0.4 

40.0 ±0.7 

86.9 ±1.0 
27.0 ±0.7 

7.3 ±0.6 
0.7 ±0.5 



0.020 ±0.002 0.020 ±0.002 4.6 ±0.5 

0.012 ±0.003 0.012 ±0.003 2.9 ±0.6 

1.000 ±0.011 1.000 ±0.011 250.1 ±1.9 

0.010 ±0.002 0.010 ±0.002 2.6 ±0.6 

1.436 ±0.014 1.435 ±0.014 381.3 ±2.3 

4.218 ±0.034 4.215 ±0.034 1154.2 ±3.9 



0.097 ±0.003 
0.007 ±0.002 
0.010 ±0.002 



0.097 ±0.003 
0.007 ±0.002 
0.010 ±0.002 



2.761 
0.016 
0.027 
0.038 
0.026 
0.032 
0.029 



±0.024 
±0.002 
±0.002 
±0.002 
±0.002 
±0.002 
±0.002 



2.733 
0.016 
0.027 
0.038 
0.026 
0.032 
0.029 



±0.025 
±0.002 
±0.002 
±0.002 
±0.002 
±0.002 
±0.002 



39.2 ±1.1 
3.4 ±0.7 
4.7 ±0.9 

1438.9 ±6.0 
8.6 ±1.1 

14.3 ±1.0 
20.9 ±1.3 

14.7 ±1.1 

19.8 ±1.2 

17.9 ±1.2 



0.005 ±0.011 
1.22 ±0.01 
0.1 ±0.3 



a in units 10 14 erg s 'cm 2 . 



are taken to be 1% (Oke |1990] , Bohlin |l99q ). These errors 
have been propagated in the calculations of the element 
abundance errors. The observed (F(X)) and extinction- 
corrected (/(A)) emission line fluxes relative to the H/3 
emission line fluxes, the equivalent widths EW of the 
emission lines, the extinction coefficient C(H/3), the ob- 
served flux of the H/3 emission line, and the equivalent 
width of the hydrogen absorption lines for region a, de- 
rived from Keck II and MMT spectra with the slit position 



P.A. = -41° are listed in Table | 
coefficient C(H/3) is small. 



The derived extinction 



The ionic abundances of 0+ 2 , Ne+ 2 , Ar+ 3 , He+ 
and He +2 were obtained using the electron temperature 
T e (Oni). The electron temperature T e (On) is adopted 
to calculate the + , N + , S + and Fc +2 ionic abundances, 
while the intermediate value of the electron temperature 
T e (Sin) was used to derive the ionic abundances of Ar +2 
and S+ 2 (Garnett [l99l ). 
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Table 3. Physical conditions and element abundances in 
the brightest H II region of SBS 0940+544 



4000 5000 6000 

Wavelength (A) 



7000 



Fig. 5. Spectra of the brightest H n region a in SBS 
0940+544 obtained with: (a) the Keck II telescope and 
(b) the MMT. The lower spectra are the observed spectra 
downscaled by a factor of 100. 



The total heavy element abundances were obtained 
using ionization correction factors (ICF) following Izotov 
et al. Q1994 |1997d) and Thuan et al. (fL99§). For oxygen 



the total abundance is the sum of 0+, CF^and O 
abundances. Although emission lines of +3 are absent in 
the optical spectrum, this ion resides in the He +2 region. 
Therefore, we derive the +3 ion abundance using the He 
II A4686 A em ission line intensity as described by Izotov & 
Thuan ( 1999 ). The ionic and heavy element abundances 
for the brightest Hn region in SBS 0940+544 together 
with electron temperatures and electron number densities 
are given in Table || along with the adopted ionization 
correction factors. 

The oxygen abundance derived from the Keck II and 
MMT spectra are 12 + log(0/H) = 7.50 ± 0.01 and 7.46 



± 0.02 respectively. For comparison, Izotov et al. ( 1991 ) 
derive d 12 + log(0/H) = 7.52 ± 0.12, while Izotov et al. 
(|!994p obtained 12 + log(0/H) = 7.37 ± 0.02. The former 
value was calculated using a three-level atom model. If the 
five-level atom model is adopted instead then the oxygen 
abundance is decreased by ~ 0.04 dex (Izotov & Thuan 



Value 



Keck II 



MMT 



T e (0 ra)(K) 
T e (0 11) (K) 
T e (S m)(K) 
N e (S ll)(cm- 3 ) 
iV e (He 11) (cm" 3 ) 
r(A3889) 

O+/H+(xl0 5 ) 
O +2 /H+(xl0 5 ) 
O +3 /H+(xl0 6 ) 
O/H(xl0 5 ) 
12 + log(0/H) 

N+/H+(xl0 7 ) 

ICF(N) a 

log(N/0) 

Ne +2 /H+(xl0 5 ) 

ICF(Ne) a 

log(Ne/0) 

S+/H+(xl0 7 ) 
S +2 /H+(xl0 7 ) 
ICF(S) a 
log(S/0) 

Ar+ 2 /H+(xl0 7 ) 
Ar +3 /H+(xl0 7 ) 
ICF(Ar)" 1 
log(Ar/0) 

Fe +2 /H+(xl0 7 ) 
ICF(Fe) a 
log(Fe/0) 
[O/Fe] 

He+/H+(A4471) 

He+/H+(A5876) 

He+/H+(A6678) 

He + /H+(mean) 

He +2 /H+(A4686) 

He/H 

Y 



18830+220 
15320+170 
17330+190 
170+60 
10+10 
0.0 

0.389+0.013 
2.727+0.081 
0.201+0.033 
3.137+0.082 
7.497+0.011 

0.874+0.034 
8.05 
-1.649+0.021 

0.476+0.015 
1.15 
-0.758+0.018 

0.608+0.015 
3.643+0.213 
2.10 
-1.546+0.025 

0.744+0.021 
1.124+0.069 
1.01 
-2.219+0.020 

0.934+0.212 
10.1 
-1.537+0.099 
0.103+0.099 



0.0801+0. 
0.0821+0. 
0.0813+0. 
0.0815+0. 
0.0005+0. 
0.0820+0. 
0.2468+0. 



0023 
0015 
0023 
0011 
0001 
0011 
0034 



20000+370 
15690+260 
18300+300 
10+10 
60+10 
0.5+0.3 

0.447+0.020 
2.438+0.102 

2.885+0.116 
7.460+0.018 

1.115+0.112 
6.46 
-1.603+0.048 

0.436+0.020 
1.18 
-0.748+0.028 

0.583+0.031 
2.794+0.516 
1.83 
-1.670+0.069 

0.751+0.050 
0.959+0.205 
1.02 
-2.219+0.059 

0.547+0.384 
8.07 
-1.815+0.313 
0.395+0.313 

0.0829+0.0066 
0.0805+0.0024 
0.0804+0.0055 
0.0807+0.0021 

0.0807+0.0021 
0.2439+0.0066 



TCF is the ionization correction factor. 



1999 ). The latter value is significantly l ower t han the pre- 
vious ones. H oweve r, Izotov & Thuan ( 1998b| ) noted that 
Izotov et al. ( 1994 ) used an erroneously low [O III] A5007 
emission line intensity. They rederived the oxygen abun- 
dance using the data of Izotov et al. (1994) and found 12 
+ log(0/H) = 7.43 ± 0.01. Thus, all 12 + log (O/H) val- 
ues are confined in the narrow range 7.43 - 7.50, or Zq/31 
- ^0/26.0 

The ratios of other heavy element abundances to the 
oxygen abundance derived from the Keck II and MMT 



12 + log(0/H) Q = 8.92 (Anders & Grevesse |TJI§) 
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observations (Table ^|) are mutually consistent and they 
are also in agreement with the heavy element abundance 
ratios derived earlier in SBS 0940+544 and in other low- 



metallicity BCDs (Izotov & Thuan |l999|) . 

The very low metallicity of SBS 0940+544 and the 
presence of very strong emission lines in its spectrum 
makes it one of the best galaxies for the determination 
of the primordial helium abundance. For this, we use the 
five strongest He i A3889, A4471, A5876, A6678 and A7065 
emission lines in both Keck II and MMT spectra. The 
first and last He I emission lines are more sensitive to col- 
lisional and fluorescent enhancement mechanisms. They 
are used to correct other He I emission lines for these ef- 



fects (Izotov et al. 1994, 1997c). The helium abundance is 
derived from the corrected intensities of the He I A4471, 
A5876, A6678 emission lines and is shown in Table 0. The 
mean values of the 4 He mass fraction Y = 0.247 ± 0.003 
(Keck II) and Y = 0.244 ± 0.006 (MMT) (see Table |) are 
consistent with the previously derived value Y = 0.247 ± 



0.007 (Izotov & Thuan 1998b) and they are close to the 
primordial 4 He mass fraction Y p = 0.244 ± 0.002 derived 
by extrapolating the Y vs O/H linear regression to O/H 
= (Izotov & Thuan |l998b| ), or to Y p = 0.245 ± 0.002 



derived in the two most metal-deficient BCDs I Zw 18 and 



SBS 0335-052 (Izotov et al. |1999|) 

It was pointed out in Sect. 2/2 that atmospheric dis- 
persion may introduce some uncertainties in the abun- 
dance determination. To estimate these uncertainties we 
measure the fluxes of the [O n] A3727 and [O ill] A4363 
emission lines in an aperture l"xl" placed on the max- 
imum in the flux distribution, and those in an aperture 
with the same size but displaced along the slit by 0'.'4 and 
O'.'l respectively for the [O n] and [O in] lines. We obtain 
relative flux differences of ~ 7% and ~ 1% for the [O n] 
and [O ill] emission lines. The effect of the atmospheric 
dispersion results in an uncertainty not exceeding 1% of 
the total oxygen abundance, because oxygen in the H n 
region is mainly +2 (see Table ||). However, because the 
H II region is much more extended than the size of the 
aperture used for analysis of the atmospheric dispersion, 
and of the large apertures used in the extraction of spec- 
tra for the abundance analysis, the effect of atmospheric 
dispersion is expected to be smaller for the [O n] A3727 
emission line and negligible compared to the statistical 
errors for the [O in] A4363 emission line. Therefore, we 
decided not to include it in the error budget. 

5. Age of the underlying stellar population 

One of the key problems debated in recent years is the evo- 
lutionary status of very metal-deficient BCDs: are these 
galaxies young or old? To resolve this problem an analysis 
of the stellar populations is required. The metallicity of 
the galaxy is one of the key parameters that determines 
the photometric properties of stellar populations. Izotov 
& Thuan ( 1999 ), considering the abundance ratios of the 
heavy elements in BCDs, have suggested that the BCDs 
with 12 + log (O/H) 7.6 might be young unevolved 




Z = 0.0004 (instantaneous burst) 

— — Z=0.001 (instantaneous burst) 

Z=0.001 (cont. SF with t f =0 Myr) 

Z=0.001 (cont. SF with t,=8 Myr) 



log t, (yr) 

Fig. 6. Equivalent widths of the nebular emission lines Ha 
and H/3 vs. age for an instantaneous burst and a heavy el- 
ement mass fraction Z = 1/20 Zq (thick dashed lines) 
and Z = 1/50 Zq (thick solid lines) calcul ated w ith the 
PEGASE.2 code (Fioc & Rocca-Volmerange |l997|) . Model 
predictions are also shown for the case of constant con- 
tinuous star formation starting at an age defined by the 
abscissa t; and stopping at ti, with if = Myr (dotted 
lines), and tf = 8 Myr (dash-dotted lines). Thin solid and 
dotted horizontal lines show the upper and lower limits 
of observed EW(Ha) and EW(Hf3) for 4 regions of the 
underlying galaxy. The shaded rectangular region shows 
the age limits in the case of an instantaneous burst. 



galaxies. SBS 0940+544 fulfills this abundance condition. 
Furthermore, the nitrogen-to-oxygen abundance ratios 
log(N/0) = -1.65 ± 0.02 and -1.60 ± 0.05 derived respec- 
tively from the Keck II and MMT spectra (Table ||) are 
confine d to the narrow range obtain ed ea rlier by Thuan 
et al. (|1995|) and Izotov & Thuan ( |l999[) for extremely 
metal-poor BCDs. The constancy of the N/O abundance 
ratio in these BCDs with a very small spread of values 
around the mean suggests that the production of primary 
nitrogen occurs in massive stars only and hence these sys- 
tems are likely to be young, since intermediate-mass stars 
have not had time to release their nucleosynthetic prod- 
ucts. On the other hand, some recent chemical evolution 
models have been succesful in reproducing the observed 
constancy of the N/O ratio, with nitrogen produced by 



the longer-lived intermediate-mass stars (Pilyugin 1999 
Henry et al. |200C| ) 



Deep V, R, I, and Ha imaging, high signal-to-noise ra- 
tio spectra showing strong emission lines and the detection 
of hydrogen Balmer absorption lines in the main body of 
SBS 0940+544 enable the study of stellar populations and 
constraint of their ages through various techniques. 

The light of the brightest H n region of the galaxy is 
dominated by a very young stellar population. The very 
blue continuum and large EW(Rf3) = 241 - 250 A (Fig. 
|], Table ||) are consistent with an age of an instantaneous 
burst not exceeding 4 Myr, for a metallicity Z = 1/20 - 
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1/50 Zq (Fig. ||). However, despite that young age, no 
Wolf-Rayet stars are detected in our high signal-to-noise 
ratio spectra. We also note that the nebular He n A4686 
emission line is not detected in the MMT spectrum and 
its intensity in the Keck II spectrum is only 0.6% of H/3, 
or several times weaker than in the very metal-deficient 
BCDs I Zw 18, SBS 0335-052 and Tol 1214-277 with de- 



tected Wolf- Rayet stellar popul ation s (Izotov et al. 1997a 
Izotov et al 



1999 



2001) 



Fricke et al 

The interpretation of the photometric data for the 
brightest H II region is not straightforward because of the 
strong contamination of stellar radiation by the emission 
of the ionized gas. The colour V — I ~ -0.5 mag is too 
blue to be explained even by the youngest stellar pop- 
ulation, because ~ 40% of the V light is contributed by 
the very strong oxygen emission lines [O in] A4959, A5007. 
Therefore broad-band photometry alone is not sufficient to 
constrain the age of young star-forming regions. However, 
spectroscopic data in conjunction with photometric data 
allows us to disentangle stellar from gaseous emission. 

Data for bright H n regions do not permit the detection 
of the small fraction of light contributed by an old stellar 
population ( e.g., P apaderos et al. 199S; Fricke et al. 2001; 
Noeske et al. 200C). The old population in such regions, if 
present, is hidden by the strong emission from the young 
stars and the gas ionized by the numerous massive stars in 
the star-forming region. More promising for the detection 
of an old stellar population is the study of the underlying 
extended emission from the host galaxy as gaseous con- 
tamination is minimized there. Nearly all BCDs show an 
underlying stellar component with red colours, consistent 
with ages greate r than a few G yr (Loose & Thuan 1986 ; 
Papaderos et al. |l996a| |l996b| Telles & Terlevich |l997| ). 
However, the colours of the extended emission in the few 
very low metallicity BCDs are consistent with younger 
ages. 

Another method for deriving ages is by fitting the ob- 
served galaxy spectral energy distribution (SED) with the- 
oretical SEDs for various stellar population ages and star 
formation histories. This method has been applied to sev- 
eral extremely metal-deficient BCDs with Z = (1/20 - 
1/40)Z Q (e.g., SBS 0335 -052 (P apaderos et al. |l99S| ), SBS 
1415 +437 (Thuan et al. |1999aj ) and Tol 1214-277 (Fricke 
et al. [2001 )). The SED of the underlying stellar component 
in these galaxies, obtained after subtraction of the ionized 
gas emission from the observed spectrum, could be fitted 
by a stellar population not older than a few hundred Myr. 

However, both methods are subject to uncertainties 
introduced by extinction. Therefore, other less extinction- 
dependent methods are desirable to better constrain the 
age of stellar populations. We use in this Section two such 
methods. One relies on the Balmer nebular emission line 
equivalent widths and the other on the Balmer stellar ab- 
sorption line equivalent widths. Such an analysis is feasible 
because of the exceptionally high-quality of the Keck II 
and MMT spectra. Because age estimates from all above 
mentioned methods depend on the adopted star formation 
history in the galaxy, we consider next different scenar- 



ios of star formation (instantaneous burst, continuous star 
formation) with varying extinction and star formation rate 
(SFR). We then put together the constraints imposed on 
the stellar age and extinction by all four methods to draw 
a consistent picture for the populations in SBS 0940+544. 

5.1. The case of an instantaneous burst 

5.1.1. Age from the nebular emission lines 

Under the assumption of an ionization-bounded H II re- 
gion, the strongest hydrogen recombination emission lines 
Ha and H/3 provide an estimate for the age of the young 
stellar population when late O and early B stars are still 
present. The ionizing flux from such a young stellar pop- 
ulation and hence the equivalent widths of the Balmer 
emission lines have a very strong temporal evolution. 

Fluxes and equivalent widths of the Ha and H/3 emis- 
sion lines are obtained for all four regions in the main 
body. Because the H/3 emission line is narrower than the 
absorption line in these regions and does not fill the ab- 
sorption component, its flux was measured using the con- 
tinuum level at the bottom of the absorption line. This 
level has been chosen by visually interpolating from the 
absorption line wings to the center of the line. The ex- 
tinction coefficient C(H/3) in those regions is derived from 
the Ha/H/3 flux ratio. For this, the theoretical recombina- 
tion Ha/H/3 flux ratio of 2.8 is adopted, which is typical 
for hot low-metallicity H II regions. No correction for the 
absorption line equivalent width has been made in this 
case. Results of the measurements together with errors 
are shown in Table ^. In Fig. ^| we compare the measured 
Ha and H/3 emission line equivalent widths with those 
predicted at a given age of an instantaneous burst. The 
time evolution of the Ha and H/3 emission line equiva- 
lent widths for the heavy element mass fractions Zq/50 ( 
thick solid lines) and Zq/ 20 (thick dashed lines) is calcu- 
lated using the gal actic e volution code PEGASE.2 (Fioc & 
Rocca-Volmerange 1997 ). This co de is based on the Padua 
isochrones (Bertelli et al. 1994) and stellar atmosphere 
models from Lejeune et al. (1998). Thin solid and dot- 
ted horizontal lines indicate respectively the highest and 
lowest observed equivalent widths of Ha and H/3 for the 4 
regions of the underlying galaxy. The shaded region shows 
that the ages range between 10 and 20 Myr, meaning that 
the gas in the main body is likely ionized by early B stars 
rather than by O stars. 

This age estimate is valid if the population of ionizing 
stars is large enough and can be described by an initial 
mass function (IMF) with an analytical form, for example 
a Salpeter IMF. However, in the case of the main body 
of SBS 0940+544, the number of ionizing stars is small 
and stochastic effects in the distribution of stars of differ- 
ent masses might be important. In particular, because of 
stochastic star formation, the ionization may be caused by 
a few O stars instead of a population of B stars. Taking the 
number of ionizing Lyc photons to be respectively 10 49 s _1 
and 10 48 s _1 for a 07v and a BOv star (Vacca, Garmany 
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& Shull 1996), we find from the observed fluxes of the H/3 



emission line that the number of 07v and BOv stars is in 
the range 2-7 and 23 - 70, respectively, depending on 
the location in the main body (see Table ^) . Assuming a 
Salpeter IMF with a slope 2.35, and upper and lower star 
mass limits of 120 M© and 2 Mq, we derive a total mass 
10 4 - 10 5 Mq for the y oung stellar population. Cervino, 
Luridiana & Castander ( 200C| ) have analyzed the stochas- 
tic nature of the IMF in young stellar clusters with solar 
metallicity. In particular, they find that in the range of 
equivalent widths EW(H(3) = 1 - 10 A, the dispersion of 
age at fixed EW(K0) is not greater than 5% - 10% if the 
total mass of the cluster lies in the range 10 4 - 10 5 Mq. 
Hence, we argue that stochastic effects in SBS 0940+544 
do not significantly change our age estimate from the hy- 
drogen emission line equivalent widths. 
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Fig. 7. Keck II spectra of regions 1 and 3 (respectively at 
6'.'3 and 12'.'2 from the brightest H n region) showing H7, 
US and other hydrogen and Ca 11 absorption lines. 



5.1.2. Age from the hydrogen stellar absorption lines 

Another extinction-insensitive method for determining the 
age of stellar populations is based on equivalent widths of 
absorption features. It can be used for larger ages than 
the nebular emission line method because the longer-lived 
B and A stars show more prominent hydrogen absorption 
lines. 



Gonzalez Delgado & Leither er (1999 ) and Gonzalez 
Delgado, Leitherer & Heckman ( 1999 ) calculated a grid 
of synthetic profiles of stellar hydrogen Balmer absorption 
lines for effective temperatures and gravities characteristic 
of starburst galaxies. They developed evolutionary stellar 
population synthesis models, synthesizing the profiles of 
the hydrogen Balmer absorption lines from H/3 to H13 for 
instantaneous bursts with ages ranging from 10 6 to 10 9 
yr, in the case of a stellar initial mass function with a 
Salpeter slope and mass cutoffs Mi ow = 1 M Q and M up 
= 80 Mq. Their models predict a steady increase of the 
equivalent widths with age. However, at larger ages 1 
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Fig. 8. The equivalent widths of H7 (a) and H<5 (b) ab- 
sorption lines in 4 regions along the main body of SBS 
0940+544. Filled large circles, open large rhombs and solid 
lines show the corrected EW & ba(Rl) and £W a bs(H<5) ob- 
tained respectively from the Keck II and MMT observa- 
tions with position angle of the slit P. A. = -41°. Large 
open circles connected by solid lines show the data from 
the MMT spectrum with the slit oriented at P. A. = 0°. 
Small symbols and dashed lines denote the equivalent 
widths uncorrected for emission. 



Gyr the situation is the opposite and the equivalent widths 
of the absorption lines decrease with age (Bica & Alloin 
1986). 

The hydrogen absorption lines due to the underlying 
stellar populations are seen clearly along the whole elon- 
gated main body of SBS 0940+544. We show in Fig. | the 
H7, H<5 and other hydrogen absorption lines in two regions 
along the slit oriented with position angle P. A. = -41°. A 
weak K Ca 11 absorption line is also present. While the H7 
nebular emission line is superimposed on the absorption 
line in region 1, it is not seen in region 3. 

We measured the parameters of the H7 and H<5 ab- 
sorption lines in 4 regions along the main body. A careful 
placement of the continuum is very important for deriving 
accurate EWs. To define the continuum level, we select 
several points in the spectral regions free of nebular and 
stellar lines. The continuum is then fitted by cubic splines 
and the quality of the continuum fit is visually checked. 
Although the contamination of the absorption lines by 
nebular emission in regions 2, 3 and 4 is very small, not 
more than a few percent, we have corrected the absorp- 
tion equivalent widths for it. For this purpose the expected 
fluxes of the H7 and US emission lines have been calcu- 
lated from the flux of the H/3 emission line. Then these 
fluxes have been subtracted from the negative fluxes of 
the absorption lines. Here we adopt theoretical recombi- 
nation Hy/H/3 and US /H/3 flux ratios at an electron tern- 



The evolutionary status of the BCD SBS 0940+544 13 



Table 4. Fluxes, equivalent widths of the Ha and H/3 emission lines and the extinction coefficient C(H/3) in different 
regions of the SBS 0940+544 main body 



Telescope 


# of region 


Distance" 


Aperture 11 


F{Ha) c 


F{HP) C 


EW(Ra) d 


EW(H/3) d 


C(H#) 


Keck IP 


1 


6.3 


1.0x3.0 


13.30 ±0.14 


3.91 ±0.12 


67.18 ±0.75 


12.37 ±0.24 


0.13 ±0.03 




2 


9.1 


1.0x2.6 


3.15 ±0.08 


0.72 ±0.08 


22.23 ±0.56 


3.26 ±0.24 


0.29 ±0.05 




3 


12.2 


1.0x3.6 


2.06 ±0.10 


0.43 ±0.08 


13.84 ±0.64 


2.18 ±0.31 


0.35 ±0.07 




4 


16.4 


1.0x4.8 


3.60 ±0.13 


1.05 ±0.12 


45.72 ±1.71 


9.15 ±0.73 


0.13 ±0.05 


MMT e 


la 


6.3 


1.5x3.0 


14.91 ±0.39 


4.22 ±0.29 


88.14 ±2.32 


14.82 ±1.01 


0.15 ±0.04 




2a 


9.0 


1.5x2.4 


3.53 ±0.17 


0.79 ±0.14 


28.76 ±1.41 


3.87 ±0.69 


0.30 ±0.06 




3a 


12.6 


1.5x4.8 


3.05 ±0.31 


0.79 ±0.28 


16.10 ±1.62 


2.82 ±0.73 


0.21 ±0.08 




4a 


16.8 


1.5x3.6 


2.96 ±0.24 


0.77 ±0.19 


41.86 ±3.41 


6.71 ±1.13 


0.21 ±0.09 


MMT' 


5 


8.0 


1.5x10.8 


14.64 ±0.57 


3.65 ±0.39 


32.74 ±1.28 


6.40 ±0.45 


0.23 ±0.06 




6 


7.0 


1.5x7.2 


14.32 ±0.47 


4.00 ±0.50 


39.75 ±1.32 


5.97 ±0.48 


0.16 ±0.07 



"distance in arcsec from the brightest H II region. 

'aperture xxy where x is the slit width and y the size along the slit in arcsec. 
c in units 10 -16 erg s _1 cm~ 2 . 
d in A. 

e slit orientation with position angle P. A. = -41°. 
■^slit orientation with position angle P. A. = 0°. 



Table 5. Uncorrected and corrected equivalent widths of the H7 and US absorption lines in the main body of SBS 
0940+544 



Telescope 


# of region 


Distance" 


Aperture 6 


EW(m) c 


EW (H 7 ) c 


EW{BS) COI c 


£W(H 7 ) cor c 


Keck II d 


1 


6.3 


1.0x3.0 


4.51 ±0.22 


2.57 ±0.23 


6.41 ±0.49 


6.37 ±0.55 




2 


9.1 


1.0x2.6 


6.96 ±0.24 


5.47 ±0.48 


7.47 ±0.68 


6.50 ±0.77 




3 


12.2 


1.0x3.6 


6.65 ±0.29 


5.40 ±0.24 


6.95 ±0.47 


6.01 ±0.29 




4 


16.4 


1.0x4.8 


5.71 ±0.53 


5.40 ±0.41 


7.19 ±0.79 


6.40 ±0.41 


MMT" 


la 


6.3 


1.5x3.0 


3.95 ±0.50 


0.70 ±0.60 


5.99 ±0.78 


5.08 ±0.68 




2a 


9.0 


1.5x2.4 


6.28 ±0.43 


4.50 ±0.47 


6.80 ±0.59 


5.61 ±0.55 




3a 


12.6 


1.5x4.8 


7.63 ±0.63 


5.23 ±0.70 


8.00 ±0.69 


6.00 ±0.74 




4a 


16.8 


1.5x3.6 


5.76 ±0.87 


5.01 ±0.52 


6.73 ±0.92 


5.81 ±0.68 


MMT e 


5 


8.0 


1.5x10.8 


5.85 ±0.53 


2.52 ±0.53 


6.86 ±0.65 


4.55 ±0.76 




6 


7.0 


1.5x7.2 


6.40 ±0.49 


2.71 ±0.52 


7.49 ±0.63 


4.92 ±0.88 



"distance in arcsec from the brightest H 11 region. 

'aperture xxy where x is the slit width and y the size along the slit in arcsec. 
c in A. 

d slit orientation with position angle P. A. = -41°. 
e slit orientation with position angle P. A. = 0°. 



perature T e = 10 4 K and an electron number density N e 
= 10 2 cm" 3 . We do not correct for extinction. Extinction 
correction would result in lower EW a bs and smaller ages. 
Hence the ages derived in this section are upper limits. 

We measured equivalent widths for only the H7 and H<5 
absorption lines. Although higher-order hydrogen Balmer 
absorption lines are seen in the spectrum of the main body, 
they are not suitable for age determination because of (a) 
the relatively low signal-to-noise ratio of the spectra at 
short wavelengths and uncertainties in the placement of 
the continuum in the blue region due to many blended 
absorption features, and (b) the weak dependence of thei r 
equivalent widths on age (Gonzalez Delgado et al. 1999[ ). 
Additionally, the He absorption line is blended with the H 
Ca 11 absorption line. 

The measured and corrected equivalent widths of H7 
and US absorption lines in the Keck II and MMT spectra 



are listed in Table |[ The errors include uncertainties in 
the placement of the continuum. In Fig. || the measured 
(dashed lines) and corrected (solid lines) equivalent widths 
of the H7 and US absorption lines are shown for the 4 
regions along the main body of SBS 0940+544. Note that 
the corrected equivalent widths of H7 and US do not show 
statistically significant spatial variations. 

Fig. |9] shows the predicted behaviour of the equiva- 
lent widths of H7 (thick solid line) and US (thick dashed 
line) absorption lines with age for an instantaneous burst 
at a mctallicity Z = 1/20 Z Q from Gonzalez Delg ado ct 
al. (|1999| ) for ages < 1 Gyr and Bica & Alloin fll986| ) 
for ages > 1 Gyr. The upper and lower limits of the cor- 
rected EW(Hrf) and EW(HS) in the 4 regions along the 
main body of the galaxy are shown respectively by thin 
solid and dotted horizontal lines. The shaded rectangular 
region gives age limits, ranging between 15 and 50 Myr. 
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log t, (yr) 

Fig. 9. Equivalent widths of H7 and US absorption lines 
vs age for an instantaneous burst with mctallicity Z = 
1/20 Z Q from Gonzalez Delgad o et al. ( |1999| ) (for ages < 1 
Gyr) and Bica & Alloin ( |1986[ ) (for ages > 1 Gyr), shown 
by thick solid and dashed lines respectively. The upper 
and lower limits of measured EW(H'y) and EW(H5) in 4 
regions along the main body of the galaxy arc shown by 
thin solid and dotted horizontal lines. Model predictions 
are also shown for the case of constant continuous star 
formation starting at an age defined by the abscissa t{ 
and stopping at if, with tf = Myr (dotted lines), and 
tf — 8 Myr (dash-dotted lines). The shaded rectangular 
region gives the age limits in the case of an instantaneous 
burst. 



Though slightly larger, this age estimate is consistent with 
that obtained from the nebular emission line analysis, par- 
ticularly since it constitutes an upper limit. Instantaneous 
burst models also predict low EW(Hj) and EW(HS) at 
age ~ 10 Gyr. However, this age is inconsistent with the 
one derived from the emission line equivalent widths. 

5.1.3. Age from the spectral energy distribution 

Useful constraints on the ages of stellar populations can 
also be obtained from the spectral energy distribution 
(SED). As already pointed out, this method is subject 
to the effects of interstellar extinction. However, the com- 
bination of the spectral energy distribution method with 
the methods discussed in Sect. |5.1.l| and 5.1.2 can be used 
for a simultaneous derivation of the age and extinction. 

To fit the observed SEDs, we use the galact ic evo lu- 
tion code PEGASE.2 (Fioc & Rocca-Volmerange |l997| ) to 
produce a grid of theoretical SEDs for an instantaneous 
burst of star formation and ages ranging between and 
10 Gyr, and a heavy element mass fraction of Z = 1/20 
Zq. An initial mass function with a Salpeter slope (a = 
-2.35), and upper and lower mass limits of 120 M© and 
0.1 M & are adopted. Because the equivalent widths of hy- 
drogen emission lines in all 4 regions along the main body 
are small, the contribution of the ionized gas emission is 
not significant. For this reason we do not include gaseous 
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Fig. 10. Extinction coefficient C(H/3), derived from the 
flux ratios of the emission lines Ha and H/3 in 4 re- 
gions along the main body. Filled large circles, open large 
rhombs and solid lines show data respectively from the 
Keck II and MMT spectra for the slit oriented at position 
angle P. A. = -41°. Large open circles show data from the 
MMT spectra with the slit oriented at position angle P.A. 
= 0° . Small symbols and dashed lines show the same data, 
but for 3 regions in larger apertures located at 7'.'6, 10'.'9 
and 14'.'6 from the brightest H 11 region with respective 
apertures l'.'0x5'/6, T/0x6'/2 and l'/0x8'/4 in the Keck II 
spectra and l'.'5x5'/6, l'.'5x6'/2 and T/5x8'/4 in the MMT 
spectra. 



emission in the SED calculations. Hence, photometric and 
spectroscopic data give us direct information on stellar 
populations if the interstellar extinction is known. 

The extinction in each region can be estimated from 
the observed decrement of the B aimer emission lines. 
Many strong hydrogen emission lines are present in the 
brightest H 11 region a (Fig. ||) which allow a precise de- 
termination of the extinction. From the Keck II and MMT 
spectra of this region we derive very low values of C(H/3), 
0.08 ± 0.02 and 0.00 ± 0.01 respectively. The Ha and H/3 
emission lines are also present in the main body, though 
weaker. We use the fluxes of these lines to derive C(H/3) 
in different regions along the slit. The results are shown 
in Table || and in Fig. |l(| We note that, in general, the 
extinction in the main body, derived from the observed 
Balmer decrement, is significantly larger than that in the 
brightest H 11 region, with a maximum C(H/3) ~ 0.3 in 
region 3 at a distance of ~ 12" (or ~ 1.5 kpc) from the 
brightest H 11 region. 

In principle, the extinction can also be derived from 
the SEDs of regions 1-4, adopting the age of the stel- 
lar population derived from the equivalent widths of the 
hydrogen emission and absorption lines. First we assume 
C(H/3) = 0. Then the observed SEDs are redder than the 
theoretical SED of a stellar population with age t = 20 
Myr (lower spectra in Fig. [TT| a - |TT|d). One of the rea- 
sons for the difference between the observed and theoret- 
ical SEDs is that interstellar extinction in the main body 
is not negligible and it modifies the observed spectral en- 
ergy distribution. To derive C(H(3) we adopt an age of the 
stellar population in each of the regions 1-4 equal to 20 
Myr, as a representative value derived from the hydrogen 
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emission and absorption line equivalent widths. We then 
compute the extinction coefficient C(H/3) to achieve the 
best agreement between the observed SED after correction 
for interstellar extinction and the theoretical SED. The 
observed extinction-corrected SEDs are superimposed on 
the synthetic 20 Myr stellar population SED for regions 1 
to 4 in Fig. [ll| a - [Ti"| d (upper spectra) . They are labeled by 
the derived values of C(H/3). These values are consistent 
with those derived from the Ha/H/3 flux ratios (Table ||, 
Fig. nff). However, despite this consistency, it is seen from 
Fig. O, that the observed SEDs (upper spectra) at A SS 
3900A are not well fitted by the theoretical SEDs. Hence, 
we conclude that the instantaneous burst models do not 
adequately reproduce the observed properties of the SBS 
0940+544 main body. 



5.2. The case of continuous star formation 



Our estimates for the stellar population age in SBS 



0940+544 made in Sect. 5.1 are based on the assump- 
tion of an instantaneous burst of star formation. Now we 
consider how the age of the stellar population is changed 
if star formation is continuous. For this we adopt a con- 
stant star formation rate in the interval between the time 
U when star formation starts and tf when it stops. Time is 
zero now and increases to the past. It is customary to con- 
sider models with a constant star formation rate, with star 
formation continuing until now, i.e. with if — 0. However, 
we will also consider models with tf > because of the 
following reason. The number of the ionizing photons in 
regions 1-4 derived from the Ha luminosity is in the 
range iV(Lyc) = (2 - 7) x 10 49 s _1 , several times lower 
than TV(Lyc) emi tted by a single star with mass ~ 100 
M Q (Vacca et al. |1996|) . Hence, an IMF with a high M up 
does not strictly apply. There are however several circum- 
stances which may explain the lack of ionizing photons: (a) 
stochastic effects on the IMF; (b) the major part (more 
than 50%) of the ionizing photons escapes the H II region 
or is absorbed by dust; (c) the IMF upper mass limit is 
truncated at ~ 30 M ; (d) star formation in regions 1 
- 4 stopped several Myr ago, i.e. tf > 0. We have con- 
sidered all of these possibilities (except for (a)) in fitting 
the data and found that the result is not changed sig- 
nificantly. Therefore, in the subsequent discussion we use 
models with M up = 12OAf and tf > 0. 

We use the model equivalent widths of hydrogen emis- 
sion and absorption lines and SEDs for instantaneous 
burst s (Fio c & Rocca-Vol meran ge 1997 Gonzalez Delgado 
et al. 1999] ; Bica & Alloin 1986 ) to calculate the temporal 
evolution of the equivalent widths of hydrogen emission 
and absorption lines in the case of continuous constant 
star formation. The results are given in Fig. § and | The 
temporal dependences of the equivalent widths of the H/3 
and Ha emission lines (Fig. ||), and the US and H7 absorp- 
tion lines (Fig. ^) are shown for continuous star formation 
starting at time t\, as defined by the abscissa value, and 
stopping at tf = Myr (dotted line) and if = 8 Myr 
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Fig. 11. Keck II telescope spectra of 4 regions in the 
main body of SBS 0940+544 on which synthetic contin- 
uum spectral energy distributions are superposed. Lower 
spectra in (a) - (d) are synthetic SEDs of a 20 Myr in- 
stantaneous burst stellar population superposed on the 
observed spectra uncorrected for extinction (C(H/3) = 0). 
Upper spectra in (a) - (d) are synthetic SEDs of a 20 Myr 
instantaneous burst stellar population superposed on the 
spectra corrected for extinction. Each upper spectrum is 
labeled by the extinction coefficient C(H/3) which gives the 
best fit of the theoretical SED to the extinction-corrected 
observed SED. 



(dash-dotted line). The equivalent widths of the H/3 and 
Ha emission lines and of the H<5 and H7 absorption lines 
in the spectrum of the stellar population formed between 
i; and tf are represented in Fig. ^ and |^ by the values of 
EW at time £;. At a fixed EW, it is seen that the younger 
the youngest stars, the larger the time interval t[ ~ tf, i.e. 
the older the oldest stars. 



5.2.1. Continuous star formation with young stellar 
population 

Can the observed SEDs of regions 1 - 4 be fitted with only 
a young stellar population continuously formed over the 
last 100 Myr? To fit the observed SEDs and derive the 
extinction in regions 1 - 4 of the main body in the case 
of continuous star formation, we consider star formation 
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occurring between tf = 8 Myr and t^ = 100 Myr. This 
model predicts EW(RS) = 1.2k, £W(H 7 ) = 6.0 A, and 
EW(Rf3) = 3.5A, EW{Ra) = 21.2k, close to the values 
observed in regions 2 and 3 (Tables | and §). However, 
the observed emission line equivalent widths in regions 1 
and 4 are larger than the predicted values. For these two 
regions, a more appropriate model is that with t\ < 100 
Myr and/or tf < 8 Myr. We show in Fig. 12 the results 
of our fitting. As in the case of an instantaneous burst, 
we adjust the extinction coefficient C(H/3) to achieve the 
best agreement between the observed SED after correction 
for interstellar extinction and the theoretical SED. The 
observed extinction-corrected SEDs are superimposed on 
the synthetic stellar population SED for regions 1 to 4 in 
Fig. [l2| a - [l2]d. They are labeled by the derived values of 
C(R(3). These values are in agreement with those derived 
from the Ha/H/3 flux ratio. Note that the synthetic SEDs 
with continuous star formation reproduce better the ob- 
served spectra in the spectral range A < 3900A than those 
calculated on the assumption of an instantaneous burst. 
This makes continuous star formation in the main body 
of SBS 0940+544, occuring during the last 100 Myr, a 
possible interpretation. 



5.2.2. Continuous star formation including old stellar 
population 

We consider next different continuous star formation sce- 
narios in which an old stellar population is present in 
the main body of SBS 0940+544. For this, we again 
adopt C(H/3) = and consider models with constant and 
varying SFRs with ages ranging between and 10 Gyr. 
Specifically, for a varying SFR, we consider two periods 
of continuous star formation with constant but different 
SFRs, occurring in the age interval < 100 Myr and > 
100 Myr. To quantify the recent-to-past star formation 
rate ratio, we introduce the parameter b = SFR(< < 100 
Myr)/SFR(i > 100 Myr). First we consider models with 
constant star formation during the whole - 10 Gyr range, 
models with b 



13 a, - 13p we show such 



i.e. models with b = 1. In Fig. 
SEDs superimposed on the observed spectra of regions 1 - 
4. It is evident that these models do not reproduce the ob- 
served SEDs. However, by increasing the parameter b we 
can fit the observed SEDs. These theoretical SEDs are la- 
beled by b > 1 in Fig. p~3] a - |l3|d. The predicted equivalent 
widths of the emission and absorption hydrogen lines are 
also in agreement with the observed ones. If non-negligible 
extinction is present in regions 1-4 then to fit the ob- 
served SEDs, the parameter b should be further increased. 
In particular, if the extinction C(H/3) = 0.28 derived for 
region 2 from the Ha/H/3 flux ratio is adopted then the 
parameter b should be as high as ~ 100 to fit observations. 

In summary, analysis of the spectroscopic data for the 
main body of SBS 0940+544 shows that the stellar popu- 
lation can be equally well reproduced by a model of con- 
tinuous star formation during the last 100 Myr, or by a 
model in which stars are continuously formed in the period 



1.0 





region 4 



C(H|8) = 0.23 



4000 5000 6000 
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7000 



Fig. 12. Keck II telescope spectra of 4 regions along the 
main body of SBS 0940+544 on which synthetic contin- 
uum spectral energy distributions are superposed. Spectra 
in (a) - (d) are synthetic SEDs of a stellar population 
continuously formed with constant star formation rate be- 
tween 8 and 100 Myr ago. These SEDs are superposed on 
the spectra corrected for extinction. Each spectrum is la- 
beled by the extinction coefficient C(H/3) which gives the 
best fit of the theoretical SED to the extinction-corrected 
observed SED. 



0-10 Gyr. In the former case, the extinction, as derived 
from the Ha/H/3 flux ratio, should be taken into account. 
In the latter fit to the observed SED is only pos- 

sible when b >5, i.e. the star formation rate in the main 
body of SBS 0940+544 has significantly increased over the 
last 100 Myr. 

5.2.3. Age of the low-surface-brightness host 

An upper limit to the age of the stellar population in SBS 
0940+544 can in principle be derived from photometric 
and spectroscopic measurements of the faint (/iy ~ 24 — 
24.5 mag arcsec -2 ) region d (cf. Fig. |lja). 

As discussed in Sect. 3.1 the V — I and V — R colours 



in this region are only slightly redder than the average 
respective values of 0.58 mag and 0.33 mag derived for 
the outskirts of SBS 0940+544 from Fig. |, thus it may be 
considered representative for the stellar LSB component. 
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Fig. 13. Keck II telescope spectra of 4 regions in the main 
body of SBS 0940+544 on which synthetic continuum 
spectral energy distributions are superposed. Synthetic 
SEDs shown in panels (a) - (d) and labeled by b = SFR(i 
< 100 Myr) /SFR(f > 100 Myr) = 1 correspond to stel- 
lar populations formed continuously with a constant star 
formation rate between and 10 Gyr. Synthetic spectra 
labeled by b > 1 correspond to stellar populations formed 
continuously between and 10 Gyr with a star formation 
rate enhanced by a factor of b during the last 100 Myr. 
These SEDs are superposed on the spectra uncorrected 
for extinction. 



The spectrum of region d extracted from the MMT #2 
spectrum (position angle P. A. = 0°) (cf. Fig. [IJb) within 
a 4'.' 8 x l'/5 aperture is shown in Fig. [bj] and labeled 
"1". Because of its very low intensity, the spectrum is 
smoothed by a 5-point box-car. No appreciable emission 
or absorption lines are seen in the spectrum, probably be- 
cause of the low signal-to-noise ratio. Therefore, we cannot 
estimate an age of the stellar populations in region d by 
extinction-independent methods, or decide on whether the 
red colours of this region are caused by enhanced extinc- 
tion or by an intrinsically red stellar population. Spectrum 
"1" , when not corrected for extinction, is fitted by a the- 
oretical SED of a stellar population continuously formed 
with constant star formation rate between 100 Myr and 
10 Gyr ago. 




4000 5000 6000 7000 



Wavelength (X) 

Fig. 14. MMT spectrum of region d in the LSB host of 
SBS 0940+544, smoothed by a 5-point box-car and la- 
beled as "1" . The spectrum is extracted from the two- 
dimensional spectrum obtained with the slit oriented at 
position angle P.A. = 0° in an aperture 4"8 x 1'.'5. A the- 
oretical SED of a stellar population continuously formed 
between 100 Myr and 10 Gyr ago is superposed on "1". 
The spectra labeled as "2" and "3" are the same spectra 
corrected for extinction respectively with C(H/3) =0.1 
and 0.15 and superposed on a theoretical SEDs of a stel- 
lar population continuously formed between 100 Myr and 
1 Gyr ago and 20 Myr and 1 Gyr ago. The latter two spec- 
tra and SEDs are offset up by 2xl0~ 17 erg s _1 cm~ 2 A~ 1 
and 4xl0 -17 erg s~ 1 cm _2 A _1 . 

However, in the presence of extinction the range of 
ages can be increased or the upper age limit can be signif- 
icantly reduced. If a value of C(H/3) = 0.1 is assumed for 
region d then the corrected spectrum (labeled "2" in Fig. 
|l4|) is fitted by a theoretical SED of a stellar population 
continuously formed with constant star formation rate be- 
tween 20 Myr and 10 Gyr ago. However, it is unlikely to 
reduce the lower age limit to values less than 10 - 15 Myr, 
otherwise emission lines must be present in the spectrum 
of region d. If a value of C(H/3) = 0.15 is assumed then 
the corrected spectrum (labeled "3" in Fig. |TJ) is fitted by 
the theoretical SED of a stellar population continuously 
formed with a constant star formation rate between 100 
Myr and 1 Gyr ago. The poor signal-to-noise ratio of the 
spectrum of region d precludes reliable estimates of the 
age for this region. 

5.2.4. Age from the surface brightness and colour 
distributions 

As emphasized in the previous discussion, photometric 
data alone do not allow to draw definite conclusions on the 
age of the stellar populations. However, they can provide 
a consistency check of the stellar population ages derived 
from the spectra. Are the ages of the stellar populations 
derived above compatible with the broad-band colours? 
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We derived V and I surface brightness and colour dis- 
tributions for the regions covered by the spectroscopic ob- 
servations at both slit orientations and compared them 
with predictions from our population synthesis modeling. 
The results of this comparison are shown in Fig. [l5|. The 
predicted colours obtained from convolving the theoretical 
SEDs with the appropriate filter bandpasses are shown by 
different symbols. The transmission curves for the Johnson 
V and Cousins I bands are taken from Bessel ( 1990 ). The 
zero points are from Bessel, Castelli & Plez ( 1998 ). 

Since the contribution of the ionized gas emission to 
the total brightness of the brightest H II region (a) of SBS 
0940+544 is high, the theoretical SED for this region has 
been constructed in the following way. We have used the 4 
Myr stellar population SED for a heavy element mass frac- 
tion Z — Zq/20. The observed gaseous SED is then added 
to the calculated stellar SED, its contribution being de- 
termined by the ratio of the observed equivalent width of 
the H/3 emission line to the one expected for pure gaseous 
emission. To calculate the gaseous continuum SED, the 
observed H/3 flux and the electron temperature have been 
derived from the optical spectrum (Tables || and ||). The 
contribution of the free-bound, free-free and two-photon 
continuum emission has been taken into account for th e 
spectral range from to 5 /xm (Aller |1984| ; Ferland |l980| ) . 
Emission lines are superposed on the gaseous continuum 
SED with intensities derived from spectra in the spectral 
range A3700 - 7500 A. Outside this range, the intensities 
of emission lines (mainly hydrogen lines) have been calcu- 
lated from the extinction-corrected flux of H/3. 

As evident from Fig. [TB^ d the observed colour of the 
brightest H n region is very blue V — I ~ -0.6 mag and 
cannot be reproduced by a stellar population of any age. 
In particular, the V — I colour equal to -0.1 mag of a 4 Myr 
old stellar population (upper open square) is significantly 
redder than that observed. However, the synthetic colour 
of a 4 Myr stellar population together with the ionized 
gas (filled square) is very similar to the observed colour. 
We also show by the lower open square the colour of pure 
ionized gas emission. It is very blue (V — I ~ -0.9 mag) 
if the observed intensities of the emission lines (Table |J) 
and the redshift of SBS 0940+544 are taken into account. 

We already pointed out that the contribution of the 
gaseous emission is small in the main body. Therefore, we 
do not take it into account and consider the colours of 
a stellar population continuously formed with a constant 
star formation rate between 8 Myr and 100 Myr ago and 
reddened by interstellar extinction. In Fig. |l5f b we show 
by filled circles the colours of four regions with extinction 
coefficients derived from the best of the theoretical fits to 
the observed extinction-corrected spectra (spectra in Fig. 
[l^ a - |l^d). Open circles show the predicted colours of a 
stellar population continuously formed between and 10 
Gyr ago, assuming enhanced star formation during last 
100 Myr, as defined by parameter b in Fig. |l3| a - [l3| d (up- 
per spectra). The agreement between the V—I colours ob- 
tained from the photometric data and those derived from 



the spectral energy distributions in both scenarios is very 
good. 

We also compare the colour of the southern region d 
(Fig. |l5|d) with predictions. Because of the noisy spectrum 
and large photometric errors for this region a broad range 
of ages can be consistent with the observed colour. In par- 
ticular, the colour of the stellar population continuously 
formed with a constant star formation rate between 100 
Myr and 1 Gyr ago and reddened with extinction C(H/3) 
= 0.15 is shown in Fig. [l5| d by the filled circle, while the 
colour of the stellar population continuously formed with 
a constant star formation rate between 20 Myr and 10 Gyr 
ago and reddened with extinction C(H/3) = 0.1 is shown 
by the open circle. They are in good agreement with the 
observed colour for this region. 

Finally, an upper age limit for stellar populations in 
SBS 0940+544 can be estimated from the V-R and V-I 
colours of the outermost regions without spectroscopic ob- 
servations. The mean colours in the outskirts of the LSB 
component of V-R = 0.33 ± 0.04 and V-I = 0.58 ± 0.03 
(Fig. [|) are compatible with those of a stellar population 
forming continuously for the last several Gyr and assum- 
ing that star formation continues until now. However, the 
spatial distributions of the old and young stellar popu- 
lations inferred from HST colour-magnitude diagrams of 
some galaxies, are different (e.g. Schulte-Ladbeck et al. 
1998). The oldest stars are seen at largest distances, where 
younger stars are not detected. It is likely that stars in 
these extended haloes were formed during the galaxy for- 
mation era and no new stars have formed since then, or 
stars have diffused to haloes from the inner regions where 
they were born. If this is the case for the outermost re- 
gions of SBS 0940+544, then to explain their observed 
V — R and V — I colours, the upper age limit should be 
significantly reduced, because no young stellar population 
is present there. In that case, the observed colours would 
be compatible with an upper age limit of < 1 Gyr for the 
stellar LSB component. 

6. Conclusions 

We present in this paper a detailed photometric and spec- 
troscopic study of the very metal-deficient blue compact 
dwarf galaxy SBS 0940+544 (Z ~ Z /27), a good young 
galaxy candidate. Photometric V and I data have been 
obtained with the 2.1m Kitt Peak telescope, and R imag- 
ing has been done with the 1.23m Calar Alto telescope. 
Very high signal-to-noise spectra in the optical range have 
been obtained with the 10m Keck II telescope and the 
4.5m MMT at two slit orientations. We have reached the 
following conclusions: 

1. SBS 0940+544 is a nearby BCD galaxy with cometary- 
like structure, i.e. a bright off-center H n region and an 
elongated main body. The V, R and I surface bright- 
ness profiles of the galaxy's main body are very similar, 
with an exponential scale length a ~ 320 pc. The V—I 
colour of the brightest H n region is very blue ~ -0.6 
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Distance along the slit (arcsec) Distance along the slit (arcsec) 

Fig. 15. (a) The V surface brightness distribution along the slit oriented at position angle P.A. = -41°. The origin is 
at the location of the brightest H n region. The regions used for spectroscopic analysis are labeled in (a) "1" to "4". 
(b) The V — I colour distribution along the slit with position angle P.A. = -41°. Filled circles are the predicted colours 
of a stellar population continuously formed with a constant star formation rate between 8 Myr and 100 Myr ago and 
reddened with an extinction coefficient as derived for each region in Fig. [l^, open circles are the predicted colours of 
the stellar populations continuously formed between and 10 Gyr with the enhancement parameters b as derived for 
each region in Fig. |l3|. The filled square indicates the colour of a 4 Myr stellar population to which has been added 
the observed ionized gas emission in the brightest H II region. The upper open square shows the colour of a pure 4 
Myr stellar population, while the lower open square shows the colour of pure ionized gas emission, (c) The V surface 
brightness distribution along the slit oriented at position angle P.A. = 0°. The origin is at the location of the surface 
brightness maximum in the main body coincident with region c in Fig. [j]. The regions studied spectroscopically are 
labeled as "5" and "6" . The location of the reddest part of SBS 0940+544 is marked as "region <f . (d) The V — I 
colour distribution along the slit with position angle P.A. = 0°. The filled circle shows the predicted colour of a stellar 
population continuously formed with constant star formation rate and reddened with an extinction coefficient C(H/3) 
= 0.15 as derived for the reddest region in Fig. [f4| (spectrum 3) while the open circle shows the predicted colour of a 
stellar population continuously formed with constant star formation rate and reddened with an extinction coefficient 
C(H/3) = 0.1 as derived for the reddest region in Fig. |IJ (spectrum 2). 



mag, due to the combined effect of the young stellar 
population and the ionized gas emission. The colours 
of the main body are much redder, without significant 
gradients. 

The derived oxygen abundance 12 + log(0/H) = 7.46 
- 7.50 is consistent within the errors with previous 
abundance determinations. The a-element Ne/O, S/O, 
Ar/O abundance ratios are in good agreement with 
the mean ratios derived from previous studies of BCDs 
(Izotov & Thuan 1999| ) . The nitrogen-to-oxygen abun- 
dance ratio log N/O = -1.60 - -1.65 lies in the nar- 



row range obtained by Thuan et al. (1995) and Izotov 
& Thuan ( 1999| ) for the most metal-deficient BCDs. 
These abundances suggest that SBS 0940+544 is a 
good local young galaxy candidate. 
The 4 He mass fractions Y = 0.247 ± 0.003 and 0.244 
± 0.007 derived respectively from the Keck II and 
MMT observations of the brightest H II region in SBS 
0940+544 are in good agreement with previous deter- 



minations of the helium abundance in this galaxy and 
are consistent with the value of the primordial 4 He 
abundance Y v = 0.244 - 0.245 derived by Izotov & 



Thuan ( |998tj) and Izotov et al. (|l999j). 
The hydrogen Ha and H/3 lines in the main body are 
seen in emission, while higher-order Balmer lines are 
seen in absorption. Three methods are used to con- 
strain the age of the stellar populations in the main 
body of SBS 0940+544. The first one is based on the 
age dependence of the Ha and H/3 emission line equiv- 
alent widths, the second one on the age dependence 
of the hydrogen H7 and US absorption line equivalent 
widths, and the third one on the age dependence of 
the spectral energy distribution. Several star forma- 
tion histories have been considered. The single instan- 
taneous burst models do not reproduce the observed 
SEDs, suggesting that star formation in the main body 
was continuous. We find that models of continuous star 
formation with a constant star formation rate starting 
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10 Gyr ago are also excluded. However, models with 
continuous star formation during the period 0-10 
Gyr, and with a varying star formation rate are able 
to explain the observed properties. In particular, mod- 
els in which the star formation rate during the last 100 
Myr has increased several times, can reproduce both 
observed equivalent widths of the emission and absorp- 
tion hydrogen lines and SEDs. However, the observed 
spectroscopic and photometric characteristics are re- 
produced equally well by models in which stars were 
continuously formed during the last 100 Myr only, if 
an extinction correction as derived from the Ha/H/3 
flux ratio is applied to the SEDs. 
5. The age of the reddest very low- intensity southern re- 
gion in the SBS 0940+544 main body gives an up- 
per limit to the galaxy's age. Because of the region's 
faintness, we were not able to detect hydrogen emis- 
sion and/or absorption features to derive a reliable age. 
Assuming no extinction, we find that the SED can be 
fitted by a stellar population continuously formed with 
a constant star formation rate for a whole Hubble time. 
However, some extinction is likely to exist. In partic- 
ular, if an extinction coefficient C(H/3) = 0.15 is as- 
sumed, the extinction-corrected spectral energy distri- 
bution of the southern region can be well fitted by a 
model with a stellar population continuously formed 
with a constant star formation rate between 100 Myr 
and 1 Gyr ago. 

In summary, we find no compelling evidence which 
favours either a young or an old age for SBS 0940+544. 
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